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We discuss in detail a five-dimensional Minimal Supersymmetric Standard Model compactified 
on S 1 fZi extended by the effective Majorana neutrino mass operator. We study the evolution of 
neutrino masses and mixings. Masses and angles, in particular the atmospheric mixing angle #23, 
can be significantly lowered at high energies with respect to their value at low energy. 
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I. INTRODUCTION 

o 

Within the Standard Model (SM), the masses of the quarks and charged leptons are determined via Yukawa 
couplings to the Higgs boson. The origin of their structure (masses and mixing angles between flavours) has no 
5h ■ explanation in the context of the SM and it is one of the major challenges for physics beyond the SM. Neutrinos being 
massless in the SM, the experimental evidence for nonzero neutrino masses gives an important indication for physics 
beyond the SM. In any case, neutrino masses are many orders of magnitude smaller than those of quarks and charged 
' leptons. Moreover, the experimental results indicate that the lepton mixing matrix has two large mixing angles unlike 
the quark mixing matrix. All this shows that the neutrino sector plays a special role in understanding the flavour 
CNJ ' structure of the SM and its possible extensions. 

Phenomcnological implications of quantum corrections to the neutrino mass and mixing parameters have been 
investigated intensively in the literature (see e.g. the review The main reason is that large effects can provide 
interesting hints for model building. From a theoretical point of view, many models are available proposing different 
explanations for the particularities of the neutrino sector. In order to study qualitative features in a model independent 
way, an attractive and simple possibility is to stick to a low-energy effective theory formulation. This means that one 
\Q . organises the effects of additional particles and symmetries present at higher energies within a systematic low-energy 
expansion. We will assume here that the heavy states arising from physics beyond the standard model completely 
decouple at low energies. In that case, the degree of suppression of an operator in the low energy effective Lagrangian 
is characterised by its mass dimension (d). The only operator appearing at dimension d = 5 is the lepton- number 
1 I violating operator Q 
Q-i! 

D . K >. 
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where L and 4> represent the lepton and the Higgs doublet fields, respectively. M is an energy scale characteristic for 
the range of validity of the low-energy effective theory description. An operator of this type can be generated, for 
instance, by the usual seesaw mechanism In that case the scale M can be identified with the mass of the heavy 
right-handed neutrino. After spontaneous breakdown of the electroweak symmetry, the Higgs acquires a vacuum 
expectation value (vev) and the operator in Eq. ([1]) then represents a Majorana mass term for the neutrinos. 
In the context of the Minimal Supersymmetric Standard Model (MSSM), it can be written in the form : 



^(L^H^L^H^) , (2) 

where L and H^ u ' now stand for the lepton and up-type Higgs doublet chiral superfield, respectively. This dimension 
five operator provides a very efficient way to study neutrino masses and mixings. Renormalisation group equations 
for this effective operator have been derived in the context of the four-dimensional SM and MSSM in Refs. 0, IE]. 

Scenarios with compactified extra-dimensions offer many possibilities for model building. For example, there are 
new ways to generate electroweak symmetry breaking or supersymmctry breaking simply by choosing appropriate 
boundary conditions (for a review on extra-dimensions and their phenomenology, see Q). In addition, for flat extra- 
dimensions the presence of towers of excited Kaluza-Klein states induces the power-law enhancement of the gauge 
couplings, leading to a possible low-scale unification 0,[U . This effect can be applied to other couplings such as Yukawa 
couplings, too, giving an original way to generate mass hierarchies For the same reasons, extra-dimensions can 

also provide a possible explanation of the observed pattern of neutrino masses and mixings. 

The aim of this paper is to study these effects explicitly in the case of one extra-dimension within a supersymmetric 
model supplemented by the effective neutrino mass operator, Eq. @. In the following we shall consider the effects of 
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renormalisation at one loop in order to test the behaviour of the extra-dimensional model. We shall focus on a five 
dimensional TV = 1 supersymmetric model compactified on the orbifold S 1 /Z 2 as a simple test ground for the effects 
of the extra dimension. 

The experimental results will be used as a starting point for the evolution of the masses and couplings in order 
to test the evolution at higher energies and the effects induced by the presence of the extra dimension. Due to the 
power-law running of the (gauge) couplings, there are of course restrictions on the range of validity of the present 
model which consequently put limits on the present investigation. 

Instead of starting from the observed masses and mixing parameters at low energies, we could take the renormali- 
sation group equations provided here, to constrain parameters of some specific model at high energies by studying the 
evolution to low energies and comparing the predictions with data. Since we are mainly interested in the qualitative 
effect of the extra dimension we refrain from adding assumptions on the high energy behaviour, although this would 
certainly be interesting from the theoretical point of view. 

The paper is organised as follows. In section 2, after a short introduction on generic TV = 1 supersymmetric 5D 
models, we present the features of a five-dimensional MSSM compactified on the orbifold S 1 jZi and discuss its low- 
energy spectrum. The details of the Lagrangian and its Feynman rules are given in appendix 1. The third section 
is devoted to a discussion of the beta functions for the Yukawa couplings and the effective neutrino mass operator. 
Numerical results for the evolution of neutrino masses and mixings are given in section 4. In the last section we 
summarise and discuss the physical implications. 



II. 5D MSSM 



A. Five-dimensional TV = 1 supersymmetry 



The beta functions can be most elegantly derived in the superfield formalism. We will therefore begin with briefly 
discussing TV = 1 supersymmetry in a five-dimensional Minkowski space and its description in terms of 4D superfields. 
More details can be found in Refs. [13, EL E2] ■ Space-time coordinates will be denoted by (x^^y). 



1. Gauge sector 

The gauge sector will be described by a 5D TV — 1 vector supermultiplet which consists (on-shell) of a 5D vector 
field A M , a real scalar S and two gauginos, A and A'. 



S a = 



d 5 a 



1 



2kg 2 



1 W MN P 
■ — Jf *MN 



D M SD M S -i\Y M D M \ 



M , 



+ -i\T M D M X'+ (A + A )[S,A + A'] 



with Dm = 8m + iAm-, T 5 = ij 5 and F MN = -y[D M , D N ]. k normalises the trace over the generators of the gauge 
groups. Equivalently, one can rearrange these fields in terms of a TV = 2, D = 4 vector supermultiplet, £1 = (V ,x), 
where: 

• V is a J\f = 1 vector supermultiplet containing A 11 and A, 

• X is a chiral TV = 1 supermultiplet containing A' and S' = S + iA 5 . 

This follows from the decomposition of the 5D supercharge (which is a Dirac spinor) into two Majorana-type su- 
percharges which constitute a TV — 2 superalgebra in 4D. Both V and x (and their component fields) are in the 
adjoint representation of the gauge group Q. Using these supermultiplets, one can rewrite the original 5D N = 1 
supersymmetric action §S§ only in terms of M = 1 D = 4 superfields and the covariant derivative in the y direction [ll[: 



S g = 



d b xd 2 ed 2 e 



4kg 



rTr 



^{W a W a S(9 2 ) + h.c) + {e- 29V V y e 29V f 



(3) 



with W a = -\D e~ 2gV D a e 2gV . D a is the covariant derivative in the 4D TV = 1 superspace (see the textbooks fl3lll4l]) 
and Vj, = d y + X- This action can be expanded and quantised to find the Feynman rules to a given order in the gauge 
coupling g . 
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2. Matter sector and its coupling to the gauge sector 

The Af = 1 supersymmetric and 5D Lorentz invariant action describing a free chiral supermultiplet is: 

S = Jd 5 x ((d M ^i)Hd M f) - ^{iV M d M + m)V) (4) 

The two complex scalars (f> ' form a doublet under a 1 SU(2)r symmetry and ip is a S'[/(2)fl-singlet Dirac spinor. 
Together, they can form two 4D Af — 1 chiral supermultiplets, $ and $ c . Adding the couplings to the gauge sector 
we obtain for the action in terms of the 4D superfields: 

S = J d 5 xd 2 6d 2 6 ($e 29V $ + $ c e - 2 f y ¥ c + ($ C (V 5 + m)<$>5{9 2 ) + h.cfj (5) 

This way of writing the action has the disadvantage that the 5D Lorentz invariance and the underlying supersymmetry 
relating $ to <I> C and V to \ are n °t manifest, but it simplifies considerably the computation of quantum corrections. 

One remark is in order here: due to the additional SU{2)r symmetry of the 5D matter sector, Yukawa-type 
couplings between cj> and ip or trilincar couplings between <E>s are forbidden in the bulk. 

B. The orbifolding and the low-energy spectrum 

If we want to recover the MSSM at low energy, we need chiral zero modes for fermions. To realize this, we will 
compactify the fifth dimension on the orbifold S 1 jZ^- The orbifold construction is crucial in order to obtain chiral 
zero modes from a vector-like 5D theory. 

The Z2 symmetry identifies y — + —y, and reduces the physical interval to [0, ttR] where R is the radius of the 
circle (see Fig. [T]). We have two orbifold fixed points invariant under the Z2 transformation, namely y = = — y and 
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FIG. 1: The orbifold projection (S 1 /Z2). The physical interval is < y < nR. 

y = ttR = —ttR + 2nR. = —y. These fixed points, called branes (cf. Fig. [2]), break the translational invariance in 
the fifth dimension and therefore the momentum conservation along y and thus part of the 5D supersymmetry. We 
will choose the transformation properties of the fields and the interactions such that 4D Lorentz symmetry, the Zi 
orbifold symmetry and one 4D Af = 1 supersymmetry are preserved. The x-field should then be odd under the Zi 
symmetry because it appears together with a derivative d y , whereas V is even. For the two matter superfields, we 
choose $ to be even and the conjugate <& c to be odd. This is a pure convention. Note that only the even fields have 



y = nR 



FIG. 2: Branes of the S 1 /Z 2 orbifold. 



zero modes. The Fourier decomposition of the fields reads: 



V(x,y) 

x(x,y) 

<5> c {x,y) 



n>l 



>1 

$(°)(o : )+a/2^$W(x)cos(^) 

n>l 



where we normalised the massive KK states to have canonical kinetic terms. At energies well below the scale R^ 1 , 
where the massive Kaluza-Klein states decouple, only the zero modes remain in the spectrum and we assume that 
physics is described by the usual MSSM. Thus, the matter superfields (and Higgs superfields) of the MSSM will be 
identified with a $' ' superfield, and the gauge fields with a mode. From the above decomposition it becomes 
obvious that the orbifold prescription breaks the original AT = 2 supersymmetry which relates $ to $ c and V to \- 

Compactifying on Si /Z 2 we have to introduce two Higgs hypermultiplets in the 4D TV = 2 language if we want to 
have zero modes corresponding to the two Higgs superfields of the MSSM, H u and Hd- Introducing an additional Z' 2 
symmetry it is possible to obtain two zero mode Higgs superfields starting from one Higgs hypermultiplct [15]. For 
simplicity we will stay with only one Z 2 symmetry here. 



C. Flavour physics on the brane 

As stated above, Yukawa couplings in the bulk are forbidden by the 5D M = 1 supersymmetry. However, they 
can be introduced on the branes, which are 4D subspaces with reduced supersymmetry. We will write the following 
interaction terms, called brane interactions, containing Yukawa- type couplings. 

Sbrane = J d 8 zdj/%) j (^\i jk ®i$j$k + | K ,.,/., //„ /,,,//„) 6(9) + h.C.j , (6) 

where $^ represents a matter superfield, L is the lepton doublet superfield, and H u one Higgs superfield. The last 
term corresponds to the effective neutrino mass operator, with dimensionful coupling kij. 

Note that we do not write any interactions with conjugate superfields < i )C on the brane which would be allowed by 
gauge interactions. But, as <& c vanishes on the branes, these interactions simply vanish. Moreover, we could have 



5 



introduced independent interactions on the ttR brane. Below we will briefly comment on this possibility, but for our 
numerical analysis we restrict ourselves to ([5]). 



III. BETA FUNCTIONS 



In this section we will derive the beta functions for the Yukawa couplings and the coupling of the neutrino mass 
operator assuming that no other operators (generated in the evolution) affect this behaviour. We will begin with 
recalling briefly the 4D result, mainly in order to set up notations and to explain the method. 



Usual 4D result 



Due to the non-renormalisation theorem [lq ]. the beta functions for the couplings of the operators in the super- 
potential are governed by the wave function renormalisation constants Zij = 1 + SZij . These relate the bare to the 
renormalised superhelds, 



N.i, 



4° = £4 /a ** ■ (?) 



The sum runs here over all iV$ chiral superfields of the model. 

In a generic 4D super- Yang-Mills theory, the result for the wave function renormalisation constant at one loop and 
in d = 4 — e dimension reads [l3| : 

1 2 ( N * 1 N " \ 

- SZy = _- £ -X* kl X jkl - 2 ]T 9lC 2 (R^)5 l3 . (8) 

\k,l=l n=l J 

We have written the interaction as Xijk^^^'^ k ' and the sum over n runs over all gauge groups of the theory. The 
group-theoretical constants C2 (R) are defined as 

C2(R)S ab = ^(T A T A ) ab , ( 9 ) 

A 

where T A are matrix representations of the generators of the gauge group corresponding to the irreducible represen- 
tation R under which the fields <&i transform. The result at two- loops can be found in [Tn |. 
For future reference we recall here the result at one-loop for the the beta function of k 1 : 

(16tt 2 )/3 k = \6Tr(Y^Y u ) - p\ ~ 6g^j k' + [Y e T Y e *] k' + k' [y}Y £ ] (10) 

The two-loop result can be found in Ref. [18j. The expression for the Yukawa couplings can be derived analogously. 



B. 5D result 



We will now derive the results in the case of the five-dimensional MSSM discussed above. To deal with the running 
in extra-dimensional theories, intrinsically non-renormalisable, we briefly remind the point of view introduced in 
Ref. The theory is treated as a chain of effective field theories where we decouple all the excitations whose mass 
exceeds the energy fi we are interested in. Hence the greater the energy, the larger the number of states considered, 
which creates to a very good approximation a power law running of the couplings (under the condition that the energy 
scale of interest is at least about an order of magnitude larger than See appendix [C] for more details. 

Higgs superfields and gauge superfields will always propagate into the fifth dimension. Different possibilities of 
localisation for the matter superfields will be studied by taking the two limiting cases of superfields containing the 
SM fermions in the bulk or all superfields containing SM fcrmions restricted to the brane, respectively. We will begin 
with the case where all matter fields propagate in the bulk. 
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1. All matter superfields propagate in the bulk 

If all matter chiral superfields of the MSSM are allowed to propagate in the fifth dimension, we find the wave 
function renormalisation constant of a matter (chiral) superfield: 

/ N g \ / JV* \ 

- (16^ 2 ) SZl D = -8 ]T glC 2 (R^)S tJ \AR + 2tt £ \* kl \ jkl A 2 R 2 . (11) 

y n=i J y k,i=i J 

Here i? corresponds to the radius of the compactified fifth dimension and A is a cutoff parameter. We only retained 
the contributions which diverge in the limit A — > oo, see appendix [C] where we discuss the evaluation of the sums over 
KK states. The same result is obtained for the Higgs superfield. A collection of explicit expressions for the different 
wave function renormalisation constants can be found in appendix ID II 

As in the previous section, the beta functions can be directly calculated from the above expression for the wave 
function renormalisation constants. Within the model discussed in Section |TT] we obtain for the beta function of 
k = k/(ttR) 2 at one loop: 

(16tt 2 )/3 k = ^(-^.g 2 - 12 5 2 )Ai? + 247rTr(yj.y„)A 2 i? 2 

+ ( [y £ t y;] k + k [y e ty e ] ) 4ttA 2 r 2 . (12) 

The beta functions for Yukawa couplings are given by: 

(i6n 2 )p Yd = y d (3Tr(yJy d ) + Tr(r e ty e ) + 3y d t y d + y tl t y t )4^A 2 i? 2 

/14 , , 32 ,\ , 

- Y d i^-g 2 + 6g 2 + -g 2 j AR (13) 

(16tt 2 )p Yu = y„(3Tr(y„ty„) + 3yjy„ + y]y rf )47rA 2 i? 2 

- Y u g 3 2 + 6. 9 2 + f 53 2 ) Ai? (14) 
(167r 2 )/3y c = Y e (3Tv(YjY d ) + Tr(Y^Y e ) + 3Y c ?Y e )4TrA 2 R 2 

- Y e (jg 2 + 6 5 2 ) AR (15) 

Note that all beta functions contain a term quadratic in AR. This term will dominate the evolution of the Yukawa 
couplings and of k. They will evolve much more rapidly than the gauge couplings, which only contain a linear term. 
This linear term arises from the sum over the tower of KK states. Since the Yukawa interactions and the effective 
neutrino mass operator are localised on the brane, we have to sum over two towers of KK excitations giving rise to the 
quadratic term. This effect has already been noticed in Ref. [l5[ , where limitations of the model due to the quadratic 
running have been mentioned, too. For the same reason at higher orders in the loop expansion higher powers of 
AR appear which limit the validity of the present approach. The top Yukawa coupling can become non-perturbative 
before the gauge couplings and at rather low energy thus limiting considerably the range of validity of the present 
model. This will become evident from the discussion of the numerical results in Section [TVl We could have introduced 
another independent interaction on the ttR brane. This would not change the general problem since it is not possible 
to mutually compensate the quadratic terms. Thus, without allowing for Yukawa interactions in the bulk, which 
would break the supersymmetry, it is not possible to avoid this quadratic running if there are matter fields in the 
bulk. 



2. Matter superfields on the brane 



In this section we will discuss the results for the beta function in the case where all matter superfields are constrained 
to live on the 4D brane at y = 0, i.e. there are no KK excitations for the matter superfields. We thus expect that 
the quadratic evolution due to the sum over two KK towers will become milder. The renormalisation constant of a 
matter (chiral) superfield becomes: 

f N g \ 
(167r 2 )<54 = -16]T 5 2 C 2 (i?«)^ +4 £ A*^, AR (16) 

V n=l k,l=l I 
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Since the Higgs-type superfields can propagate in the bulk, in contrast to the matter superfields, the wave function 
renormalisation has no longer the same structure. For example, the contribution containing the vector superfields 
contains a sum over KK states for the Higgs superfield, whereas only the zero mode enters in the case of the matter 
superfields. For the Higgs superfield we obtain 

/ N g \ / N * \ 

- (16tt 2 ) SZ H = 1-8 ]T 9 2 n C 2 (R^)8 tl \ AR + £ \* kl \ jkl log AR (17) 

V n=l J \k,l=l J 

Again, explicit expressions for the wave function renormalisation constants can be found in appendix ID 21 

As in the two preceeding sections, we can derive the beta functions from the above equations for the wave function 
renormalisation constant. The beta function of n is given by 

(16tt 2 )(3 k = ((-y<?i 2 - 18<? 2 2 )Ai? + 6Tr(Yj.Y„)^ K + {kY$Y b + Y e T .Y» 4 Ai? . (18) 

The beta functions for Yukawa couplings read 



(16n 2 )f3 Yd = Y d (3Tr(Y}Y d ) + Ti(YjY e )) 



Y d I - %g\ - 9g 2 ~ ^g 2 + \2Y\Y d + AY^Y U ) AR ( I .!) ! 



15 J1 ^ 3 
(16^ 2 )/3y„ = 3Y u Tr(YjY„) 



+ Y u (-^g 2 - 9gt - y <?I + 12Y tYu + 4Y}Y d ^j AR (20) 
(167r 2 )/3y e = Y e (3Tr{YjY d )+Tr(Y^Y e )) 

+ Y e (-^g 2 - 9g 2 2 + 12YX^) Ai? (21) 

As can be seen from the above equations neither k nor the Yukawa couplings contain a term quadratic in AR any 
more, since there are no KK excitations for the matter superfields. This means that the range of validity of this 
model will be wider since the couplings, in particular the top Yukawa coupling, will not become non-perturbative at 
low energies (a few R^ 1 ). This will be discussed in more detail in the next section where numerical results will be 
presented. 



IV. NUMERICAL RESULTS 



Within this section we will apply the beta functions derived in the previous section to study the evolution of the 
different couplings. We are particularly interested in the evolution of the neutrino mass parameters. To that end we 
employ the MATHEMATICA package REAP [l9j in a slightly modified version including our models. The SUSY 
threshold is taken at 1 TeV. Throughout the discussion we will identify the cutoff A with the energy scale /i for the 
evolution. 

Before coming to a detailed discussion of the numerical results, let us mention that the discussion of the fixed point 
structure for the mixing angles follows closely the MSSM one [2(| Hlj . The evolution equations (given in detail in 
Appendix [A} have the same structure as in MSSM. We will come back to this point in section HVB1 



A. Matter superfields in the bulk 



Let us begin with the model described in Sec. lIIIB II The corresponding beta functions have been given in Eqs. (|12l 
I15p . In this case, as already anticipated, the evolution is dominated by the quadratic terms. This puts stringent limits 
on the model. As soon as the energy scale /j, becomes of the order of i? -1 , i.e. the couplings start to feel the effect 
of the extra dimension, y t starts to increase rapidly and it diverges for ~ 2 — 3. Moreover cubic and higher 
divergences arise at higher orders in the loop expansion which can strongly affect the results. This is due to the 
presence of sums on the KK excitations in the loops which are not restricted by any conservation of KK numbers at 
the vertices. 
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Throughout the calculation of the beta functions, we have assumed that \iR 3> 1. How reliable is this approximation 
if we reach the perturbative limit of the model already at fj,R ~ 2 — 3? We checked the validity of the approximation 
by integrating the beta functions without using the power law approximation. This can be done numerically, adding 
at every threshold the contributions from the new degrees of freedom. The results do not change considerably. For 
instance, y t diverges for a slightly larger value of ~4-5 (see Fig. [3]). 



Yt 9i 




H (GeV) 



FIG. 3: Evolution of the top Yukawa coupling (left panel) and the gauge couplings (right panel) g\, Q2 and g% (from top to 
bottom) in the case ~ 10 10 GeV and matter fields in the bulk. We clearly see that yt diverges before any perturbative 
unification is possible. 

However, a unification of the gauge couplings before the non-perturbativity of the top Yukawa coupling is possible 
at a scale R^ 1 > 2.10 14 GeV. This scale is of the same order as the standard 4D unification scale, such that the 
extra dimensional scenario looses much of its interest. In addition, the rapid divergence of the top Yukawa coupling 
prevents us from making any reliable prediction on the effects of the extra dimension on the neutrino sector. We will 
thus not show any results for the neutrino parameters within this model. 

B. Matter superfields on the brane 

If we restrict the matter fields to the brane, we no longer face this problem. First of all, there is no divergent 
quadratic term at ldoop. In addition, thanks to the large negative contribution of the gauge couplings to (3y u , yt 
decreases. This allows for a perturbative unification of the gauge sector at a value fi ~ 40i?~ 1 = Mqut for R^ 1 = 10 3 
TeV. We have also checked explicitely that 2-loop terms are at most quadratic in the cut-off. These terms are further 
suppressed by an extra factor of 167T 2 so that the results are presumably not modified within the range of validity of 
the effective theory established at ldoop. 

This model is thus much more promising and worthwhile studying the evolution of the neutrino mass parameters. We 
discuss our conventions for the mixing matrix and masses together with the explicit renormalisation group equations 
(RGE) in appendix O 

Here a few general comments are in order. The RGE for the neutrino masses and mixings are similar to the 4D 
case since the beta functions have a similar structure. Consequently the relevant parameters will be essentially the 
same as in the 4D case, tan (3 plays an important role as all the mixing angles and phases depend on y T . In addition, 
mi is important since the running of Q\2 will be stronger if m\ ~ m^. This is the reason why we generally choose 
tan/3 ~ 50 and mi ~ 0.1 eV to explore the effects of the extra-dimension. 

In Fig. [4] we show the influence of tan/3 on the evolution of 6*12, all other parameters are kept fixed. Depending on 
the value of tan /3, #12 can assume almost any value at the scale where gauge coupling unification is achieved within 
our model. It indicates the scale where new physics will come into play. This sensitivity of Q\i to tan (3 has already 
been noticed in 4D (cf for instance Ref. (22J). 
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FIG. 4: Running of i2 for the values tan/3 = {10, 30, 50} (from top to bottom), i?" 1 = 10 4 GeV, m l = 0.1 eV, 8 13 = and all 
phases vanish at Mz- 

1. Masses 

Still, the 5D case is intrinsically different from the 4D one and has some particularities that are quite independent 
of any choice of low energy parameters. This is the case for the evolution of the masses (Fig. HJ. As can be seen in 
Fig. [5] the evolution has exactly the same form for the three masses and is much sharper above /i ~ R~~ , i.e. where 
the masses start to feel the effect of the fifth dimension. This leads to a reduction of up to a factor of the order of 
five for the masses in the UV with respect to the values at low energies. This prediction is extremely stable and can 
be explained as follows. 

The evolution of the masses is governed by the following equation 

d( f =kU) = 7hi = is? [a + mi ■ (22) 

where the parameters induce a priori a non-universal behaviour and the parameter a is detailed in appendix [XJ in 
4D and 5D. In contrast to the MSSM, the evolution in our case is completely dominated by the universal part. The 
essential point is that in the MSSM the positive contribution to a, approximately 6y 2 , is of the same order as the 
negative contribution from the gauge part, leaving an 0(1) factor. In the setup we are interested in, however, the 
situation is completely different: on the one hand yt decreases very rapidly and on the other hand the contribution 
from the gauge part to a is multiplied by 12/ii? with respect to the MSSM which makes it completely dominant 
compared to any other contribution. We can therefore write: 



1 ^ 18 2 1Q 2 







^tiR [ -—gi - 18g 2 2 ) mi . (23) 



This equation is universal. From this approximation we immediately see that all masses decrease with increasing 
energy and eventually become zero. 

This discussion can be extended to Am 2 , too, except in the case where mi is not small compared with m 2 and 
tan/3 is large. In the latter case we can check from the analytic formula, 

^Am 2 oi = ^ [aAm 2 sol + ^Ry 2 T [2s 2 23 (m 2 2 c 2 12 - m\s 2 l2 ) + O(0 13 )}} (24) 



that for */ Am 2 ol <C mi the first term does not necessarily dominate. Still, because all the masses decrease rapidly, 

"sol 



Am 2 ol has to stop growing at some point. This is indeed the case, see Fig.[6l 
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FIG. 5: Running of the three masses for the values i?" 1 = 10 4 GeV, mi =0.1 eV, #i3 = and all phases vanish at Mz- 




FIG. 6: Running of Am' ol and AmJ lm for the values tan/3 = 50, R 1 = 10 4 GeV, mi = O.leV and all phases vanish at Mz- 

2. Mixing angles 

Note also that, assuming a hierarchical spectrum (mi = O.OOleV) and tan/3 = 10, #12 will vary twice more than in 
the MSSM between Mz and Mqut, provided #13 is rather large (~ 9 degrees). 

Another interesting observation is that 823 can become much smaller than in the 4D case. This allows for a 
unification of 612 and 823 at a value below ten degrees. This can be seen in Fig. [7] where we display the three mixing 
angles for tan/3 = 50 and mi = O.leV. 

In Figs. [5] and [5] we show results for 823 at the cut-off scale in our 5D model and in the 4D MSSM. Fig. [5] 
corresponds to an inverse hierarchy whereas Fig. [9] corresponds to normal hierachy. For the results shown we have 
O.OleV < mi < O.leV. For smaller masses the spread is reduced. In general, compared with the MSSM we obtain new 
interesting points where 6*23 is small, (< 20 degrees) for quite large tan/3 ~ 30 — 50 and mi ~ O.leV. Note that it is 
even easier to get such small 823 with an inverted hierarchy. In this case we can also find a portion of the parameter 
space where 6*13 can be rather large, see Fig. [TDJ This is in sharp contrast with the 4D MSSM, thus suggesting the 
possibility that #13 ~ 812 at the cut-off scale with values around roughly 30 degrees. There arc even some points with 
larger values, up to 45 degrees. 

Let us note that, although previous studies [23| have investigated the possibility of having CKM-like values for 
the PMNS angles at high energies, they required even larger values O.leV < mi < 0.6eV creating some tension with 
cosmological bounds. Here a value mi < O.leV is sufficient, which is new compared with the 4D case. This opens the 
interesting possibility to find a set of parameters which generate a maximal mixing from a small CKM-like mixing at 
some new physics scale or even allow for unification of the three mixing angles at high energy. 

Despite these interesting features we note that the region in parameter space most studied here is the region where 
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FIG. 7: Running of 812 (red/grey), #13 (blue/dark grey) and 9; 
mi = O.leV and no phases at Mz 



623xd 




(green/light grey) for the values tan/3 = 50, R 1 = 10 4 GeV, 



923mssm 




15 20 25 30 35 40 45 5Q ta 



FIG. 8: Comparison of #23 at the cut-off scale as a function of tan/3 in our 5D model and in 4D MSSM for random phases and 
0.01 < mi < 0.1 eV with inverse hierarchy. For smaller values of mi the spread is reduced. 



mi ~ TO2- In this case the low energy values of the mixing angles (in particular 812) are very sensitive to the precise 
value of mi at the high energy scale (at the percent level) . 

As mentioned above, we have in principle the same fixed points for the angles as in the 4D MSSM. In the case 
where 613 = it can be seen immediately from the equations that #23 = 9 12 = are fixed points. For nonzero #13 the 
fixed points are less obvious, but as the structure of the equations does not change with respect to the MSSM, the 
result is in principle the same. It leads to fixed points predicting the relation [2^, Hl[ 

sin 2 2gi 2 = '| 3S 2 iD ?t% with S l 3 ^Ul±^l-sm 2 262 3 ) (25) 

l s 23 c 13 + s 13) * V ' 

for different patterns for the coefficients. The coefficients change with respect to the 4D MSSM and therefore the way 
we can approach the fixed points. Two essential differences between the 5D and 4D MSSM case should be noted in 
that context. First, the coefficient C governing the evolution, see appendix [Al is C = 1 in the 4D MSSM whereas it 
is 4/iR in the 5D case. Second, the numerical coefficients entering the equations for the angles depend on the masses 
and become almost constant if we enter the universal regime for the running of the masses discussed above. This can 
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FIG. 9: Same as Fig. [8] for normal hierarchy. 
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FIG. 10: Same as Fig.lfor 6>i 3 . 



be understood as follows. Neglecting the phases, these coefficients can be written as (r/ij + nij)/{rrii — rrij) such that 
the indirect dependence on fi via the masses cancels out within the universal regime. From our numerical analysis 
we find that in general the fixed points are not reached within the range of validity of the present model. Increasing 
tan/3 accelerates the evolution such that the angles come closer to the fixed points at the cut-off scale. For very large 
values, tan/3 > 50, we reach in some cases a fixed point. 



V. CONCLUSIONS 



The aim of this paper was to perform a study of the neutrino sector in a simple supersymmetric extra-dimensional 
framework. We have presented beta functions for Yukawa couplings and k, the coupling of the effective neutrino 
mass operator, within two distinct setups in the context of a five-dimensional MSSM on S /Z2- In the first case, all 
fields associated with the SM fermions are allowed to propagate in the fifth dimension whereas in the second they are 
restricted to the brane. 

Due to the 5D TV = 1 supersymmctry, Yukawa interactions are forbidden in the bulk and must be introduced on 
the branes. Within the first model, this induces a quadratic running for Yukawa couplings and n. y t then becomes 
non-perturbative already at rather low energies and even before the gauge couplings. This strongly limits the range 
of validity of the model. 

Within the second scenario the dependence on the energy scale is only linear and Yukawa couplings remain pertur- 
bative until gauge coupling unification. The evolution of neutrino masses and mixings shows interesting possibilities 
to explain the observed values at low energies from some specific scenario at high energies. As a generic prediction, 
neutrino masses are reduced up to a factor of five at the unification scale with respect to their values at low energies. 
From a top-down point of view, one can radiatively generate a large mixing pattern at low energy starting with small 
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values of 6*23 and 6*12 at a high energy scale with values for mi consistent with cosmological bounds. It is also possible 
to generate a small 6* 13 at Mz from large values at the cut-off scale. 
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APPENDIX A: CONVENTIONS FOR NEUTRINO MASSES AND MIXING PARAMETERS 

Within this section we would like to stress our conventions for the mixing angles and phases and briefly discuss 
different scenarios for neutrino masses. The mixing matrix which relates gauge and mass eigenstates is defined to 
diagonalise the neutrino mass matrix in the basis where the charged lepton mass matrix is diagonal. It is usually 
parameterised as follows [24l |: 

U = diag(e^ , e i5 " , e a - ) U M ns , with 

/ c 12 c 13 s 12 c 13 s 13 e~ lS \ 

Umns = -si 2 c 23 - a 2 s 23 si 3 e lS ci 2 c 23 - s 12 s 23 si 3 e tS _ s 23 c 13 diag(e~ I0l/2 , e~ 4<fc/2 , 1) . 
\ S12S23 - ci 2 c 23 si 3 e lS -C12S23 - si 2 c 23 si 3 e lS c 23 a 3 J 

(Al) 

where Cy = cos dij , s y - = sin 6>y . We follow the conventions of Ref. [22| to extract mixing parameters from the MNS 
matrix. 

Experimental information on neutrino mixing parameters and masses is obtained mainly from oscillation experi- 
ments. The most simple interpretation of these oscillation data is in terms of massive neutrinos. There are essentially 
three types of experiments providing us with data: solar neutrino experiments (Kamiokande and Superkamiokande) , 
looking for a deficit of v e from the sun, atmospheric neutrino experiments looking for a deficit of and produced 
by cosmic rays in the earth's atmosphere and reactor experiments looking for the neutrino flux from a reactor. Present 
data are compatible with oscillations between the three known neutrino flavours. In general Am, jB is assigned to a 
mass difference between v 3 and v 2 whereas Am 2 so[ to a mass difference between v 2 and v\. Current values [25| are 
summarised in Table lAl Data indicate that Am^ <C Am 2 tm , but the masses themselves are not determined. Either 
they follow the hierarchical scenario with Am 2 tm s» m 2 3> m 2 3> m 2 or (partly) degenerate scenarios with 

3 2(1) 1(2) 

masses approximately equal. 



Parameter 


Value (90% CL) 


sin 2 (2(9i2 ) 


0.86(i° ;°oi) 


sin 2 (2(9 2 3 ) 


> 0.92 


sin 2 (20i 3 ) 


< 0.19 


Am 2 oi 


(8.0±g;|) x lO^eV 2 




1.9 to 3.0xl0 -3 eV 2 



With these conventions we can derive the approximate equations for the mixing angles and masses. The derivation 
is similar to the one in ref. [22] , and with the same notations we write the result : 
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1 



m i 



m 2 = 



m 3 = 



16tt 2 
1 

16tt 2 
1 



[a + Cy 2 (2s 2 2 s 2 23 + 0(6 13 ))} mi 
[a + Cy 2 (2c 2 2 s 2 23 + 0(e 13 ))}m 2 



[a + Cy 2 x 2c 2 13 c 2 23 }m 3 

•71 ~ 

2 sm(26'i2)s 2 3 



32tt 



sin(26» 12 )sin(26l 2 3)- 



m 3 



0(M 



x [mi cos(</>i — (5) — (1 + C)™2 cos(02 — 5) — C TO 3 cos <S] + O{0i 3 ) 



723 — 



32tt 2 



sin(20 23 ) 



1 



Am 2 , 



c 12 \m 2 e 



m 3 \ 



2 |mie 



m 3 



i + C 



•0(013) 



where Am 2 



m 2 2 -m\, Am 2 at 



of a and C, which are in the model with matter superfields on the branc: 



and C = Am 2 oi /Am 2 tm . The main difference here lies in the expressions 
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--g 2 -18g 2 )»R + 6Tr(YlY u ) 



C = C{p) = 4nR 



(A2) 



This has to be compared with the MSSM coefficients: 



a= [6Tr(YlY u )~-g 2 -6g. 



C= 1 



(A3) 



APPENDIX B: COMPLETE LAGRANGIAN AND FEYNMAN RULES 

Here we write the complete 5D action of the model where all fields can propagate in the bulk [l( 



Sn 



5,, 



Tr 



C 2 (G) 



<I' N :<!// { -VU{P T - i(P a + P 2 ))V - VdjV - 



X"FTX + oXX 



2^D 



y -{D 2 D a V)[V, D a V] + g (d 5 V[V, X + x] - (x + xW, X + X]) 
0(g 2 ) + ghosts X 

d s zdy + $9#= + ^d 5 $i8(§) - ^8^6(9) 

g(2$iV$i - 2<S>lv$>1 + $ c iX ®it(d) + *<X^(0))} 



Sbrane 



d ,s :d.</<M.</) { ( ~\ijk$i®j$k + -k i jL i H u L j H u ) + h.c. 



where we separated the pure gauge, the coupling of matter to gauge and the Yukawa sector. The latter is localised 
on the brane in order to respect 5D SUSY. £ is the gauge fixing parameter. 

In order to consider the 4D effective theory we compactify this action by expanding the fields as in section [IT] and 
we keep only the terms that will be of interest to renormalise the Yukawa beta functions: 
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■> 'gauge 



+ V — V {n) V {n) + - V — Y (n) -Y (ll) + - V Y (Tl) Y (n) 
n>l n>l n>l 



5 matter = f d 8 z J $(0)$(0) + ^($(")$(«) + $c(n) $ c(n)) _ ^ » ($^(")$(")«5(6») + $(")$ c (™) < 5(6»)) 

n>l n>l 



n>l 

2 ^ y(o)$c(n) _|_ ^($ c (™) Y ( n )$(°)(5(0) + $(°)^(™)$ c (")(5(6»)) 



n>l 



n>l 



+9 



,n>l 



m,n>l 



_ _L (<j)C(m)^(n)^(m+7i) _ <j>(l™-™l))<J(0) _|_ $("») ^(™) ($c(m+n) _ $c(|m-n|)^^^ 



m,n>l 



5*f, r 



*ijk 



$(0)$(0)$(0)+3^2$W$( )$(°)+6 E <lf>$< m >#< n 

n>l m,n>l 



+2V2 E $i ro) ^ n) $i p) 

m,n,p>l 



<5(6») + h.c. 



+- 



+4 £ (L^FWLfFW+L^WLfffW) 



m,n> 1 



+4V2 ]T (L^H^LfH^+L^H^LfH^) 

m,n : p> 1 



+4 ]T L^H^LfH^ 

rn,?i,p,q> 1 



<f(0) + h.c. } 



We kept the interaction terms involving only excited KK modes although they are not relevant for the one-loop beta 
functions, since they can be relevant to the renormalisation of the Kaluza-Klein masses and to the mixing of the 

excited states. We defined the 4D effective couplings: g = Jj A = — , ^ k 



Gr' (v^R) 3 ' i^R) 2 ' 



The propagators are read off from the action (we choose the gauge £ = — 1) 
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j. (c)(n) / 


-p, 9) 


P 


i(c)(m)/ /j/x 




-p, 9) 


P 






-p,9) 


P 






-p,9) 


P r ^ 
"^^^^^ 0') 




-p,9) 


p 

o o o o o^^*-o o o o 





2 t ^2 — dijdmn l 5 4: {9 — 9') 

D\p) 



D\p) 



SijS mn S 4 (9 - 9') 



R 2 

2i 



V 2 + w - *e 



P 2 + W~ ie 



5 ab S mn 5 4 (9 - 9') 
SabS mn S A (9 - 9') 



the vertices that are relevant for our one loop calculation: 

$(0) 



ya(0,n) 



2igT a x a(n) 



$(0,n) 



$(0) 

igT a 

$c(n) 



and we derive 



(0) 



(0) 




(0) 



(n) 



(0) 




(0) 



(m) 



(0) 




(n) 



We can do the same for the case where all superfields containing SM fermions are restricted to the brane. The part 
of the action involving only gauge and Higgs fields is not modified, whereas the action for the superfields containing 
the SM fermions splits into: 



Strane = J d 8 zdyS(y) + 2g^V^ } 



n>l 



a 



d 8 zdy%) |f e S c ii/ d + Y d D c QH d + Y U U C QH U + hlH u LH u + h.c. 
J d s z ^Y e E c LH^ + Y d D c QH d Q) + Y U U C QH^ + ±kLH™LHW 
+ ^ (Y e E c LH { d n) + Y d D c QH d n) + Y U U C QH^ + \kLH { ^LH^ 



n>l 



^ ^Li^LtfW+h.c.l 

n , n > 1 J 
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n > 1 




n > 1 

FIG. 11: One-loop diagram with two Kaluza-Klein states running in the loop restricted to have the same KK number 

Y K 
where we have written: K = et k = . 

The propagators and Feynman rules can be derived in the same way as before. 

APPENDIX C: KALUZA-KLEIN INTEGRALS 

We will give here the major steps to the derivation of the relevant Kaluza-Klein integrals and the computation 
of divergences following [7]. We will not show finite terms, i.e. terms which vanish in the limit A — > oo. For 
the calculation of the one-loop contributions to the wave function renormalisation constants we need four types of 
two-point functions which we will discuss in turn. 

1. Two excited KK modes with the same KK number 

The first case contains two excited KK modes, but their KK number is restricted to be the same. It is illustrated 
in Fig. [IT] and arises typically from bulk interactions. For example, it enters the contribution from the one-loop 
correction containing one vector superfield and one chiral superfield (which can be only Higgs for the model discussed 
in Sect. IIIIB 21 and Higgs or matter superfield for the model in Sec. IIII B"Tj) . 
Let us start from the following general expression 

K _ y f d'k 1 1 rrn 

I {2^ k 2 ~m 2 ~n 2 /R 2 k 2 ~ M 2 -n 2 /R 2 ' 1 ' 



n>l 



After introducing Feynman parameters and performing a Wick rotation in the standard way, the resulting momentum 
integration can be re-expressed using a proper-time regularised form. This allows to introduce two cutoff parameters, 
tin and £tjv to treat infrared and ultraviolet divergences, respectively. Assuming Mq ~ too <C R~ l the sum over KK 
states can be evaluated. We obtain 



Jo Jt 



1 rtiR 



1 „ / it \ 1 



2 3 V ttR 2 J 2 



The #3 function arising from the sum over KK states is defined as: 

n— +oo 



(C2) 



,(*)= E e " nx - ( C3 ) 



n— — oo 



As discussed in the appendix of Ref. [3] we can use the approximate form of the #3 function, 



to evaluate the integral. This form of the 9% function can be applied in general if tm, £tjv ^ but it gives also a 
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very good approximation to the integral if tm ~ R . We then obtain : 

2 fl 

/o Jtvv 
ft 



(16tt 2 )/C 



dt 



I 2t 3 / 2 2* 



With the redefinitions iuv = 2 , <ir = ri? 2 and r = ir/4 the final result reads: 

(167r 2 )/C = -2 + 2Ai?- ilogA 2 i? 2 
~ 2AR- ilogA 2 i? 2 . 



(C5) 
(C6) 

(C7) 
(C8) 



In the last line we have supposed that AR ^> 1. 

In Ref 26], the authors discussed another coherent cut-off regularisation scheme which allows for obtaining the 
preceeding result (|C8j) naturally without any rescaling of the cut-off when the KK tower is truncated at AR. 



Two KK excitations with different KK numbers 



We now perform the integral where two Kaluza-Klein states run in the loop and are not constrained to have the 
same KK number, 



E 

n , m > 1 



d 4 k 



i 



l 



(2tt) 4 k 2 - ml - n 2 /R 2 k 2 - M 2 - m 2 /R 2 



(C9) 



This type of integral arises only in the model discussed in Sec. IIII B 2| where all matter superfields propagate in 
the bulk. It appears in connection with the Yukawa interactions restricted to the brane. This type of function is 
illustrated in Fig. [12] 



n > 1 




m > 1 

FIG. 12: One-loop diagram where two Kaluza-Klein states are not constrained to be equal 



Following the same steps as in the first case, we have: 



(16^ 2 )£ 



n , m > 1 



dt 



.-A nm t 



dx 
o Jt 



dt 



1 / itx 



I ( it(l-x) 
; ! irR 2 



where 



A nm = (ml + n 2 /R 2 )x + (M 2 + m 2 /i? 2 )(l - x) ~ n 2 x/R 2 + m 2 (l - x)/R 2 



In exactly the same way as in the previous section we obtain in this case 

(16tt 2 )£ = 7rA 2 i? 2 -7r + 4-4Ai?+ilogA 2 i? 2 
~ 7rA 2 i? 2 -4Ai?+ilogA 2 i? 2 



(CIO) 
(Cll) 

(C12) 
(C13) 
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3. One KK mode running in the loop 

The third type of integral contains one zero mode and one excited KK mode: 

d 4 k 1 1 



<J (27T) 4 fc 2 



i 2 /R 2 k 2 -M 2 



(C14) 



It is illustrated in Fig. [T2J The evaluation proceeds in exactly the same way as before. We obtain 

dt 



-n 2 x/R 2 t 



t 



r tiR dt 




f itx \ 


1" 


/t uv t 


> 


VTTi? 2 ) 


_ 2 



which gives upon performing the same approximations as before 

(16tt 2 )W = 4AR-4- ilogA 2 i? 2 
~ 4Ai?- -logA 2 i? 2 . 




(C15) 

(C16) 
(C17) 







FIG. 13: One-loop diagram with only one Kaluza-Klein state in the loop. 



4. Only zero modes running in the loop 



We recall the result of a loop with only zero modes, 

d 4 k 1 



(2tt) 4 fc 2 - m 2 k 2 - M 2 ' 



(C18) 



which reads 



( \ (>7t~ iJ = / cfe / * e -*(™§«+^(i-*)) ~logA 2 Jl 2 . (C19) 



o 



APPENDIX D: RENORMALISATION CONSTANTS 



Within this section we summarise the explicit expressions for all wave function renormalisation constants needed to 
compute the beta functions. In this section we will display the diagrams entering in the one-loop renormalisation of 
the Yukawa couplings and their values, in order to deduce the wave functions renormalisation and the beta functions. 
Extensive use is made of the integrals calculated in the last section. 
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1. Matter fields propagating in the bulk 



We have 5 types of diagrams: 
(n) 



(n) 
1 2 



(0) 




(n) 



(n) 



(0) 



(0) 



(m) 



(n) 



(n) 



We provide some steps of the calculation for the first diagram and give the result for the four others : 



'2(P + ^) 



1 

16' 



-iSf 2 



x -.DiDl - :; 12 „ 2 ^ (0) (-p, ^)^ (0) (P, 02) 



= -i2g z C 2 (R)S™ AC 



d 4 fc 



1 



(2tt) 4 k 2 { P + k.y 



d 4 «: (o) (-p,«i)r(p,ei) 



(0), 



= z- 



-2ff 2 C 2 (i2)<Sr 



16tt 2 



(2A J R + log(A J R)) J tfeifi {Q \-p,6Mf>{p,e{) 



sz. 



13 
(4) 



-2. 9 2 C 2 (i?)<5.; 



7 (3) 



— (2Ai? - log(Ai?)) ^ 



. 2AifejA* /c; <5 rs 



(4AJ?-log(A J R)) 5Z, 



(5) 



171 

Ait; A*j,,(5. 



log(AiJ) 



^ikl A jkl°rs 



(27r(Ai?) 2 - 8Ai? + log(Ai?)) 



16tt 2 v fav " y 16tt 2 

We only displayed the integral over the 9 coordinate for the first contribution, omitting it in the others. 

Summing the different contributions, every subdominant divergence disappears for both the gauge and the Yukawa 
contribution and we obtain: 



N g \ / JV* 

- (16^ 2 )<5Z| D = ( -8^ 3 2 C 2 (i?«)^ ] AR+ ( 2ir ]T X* kl X jkl ] A 2 i? 2 

k,l=l 



(Dl) 



n=l 



1 
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Applying it to the matter and Higgs supcrfields : 



{l6ir 2 )5Z Hu = 12ttTt(Y^Y u )A 2 R 2 - (-g 2 + 6g 2 )AR 

5 

-(16tt 2 ) SZ Hd = 4ir{3Tv(YjY d ) + Tt(Y^Y c )}A 2 R 2 - {p 2 + 6g 2 )AR 

-(16tt 2 ) SZ L = 4ir(Y e ^Y e )A 2 R 2 - £g 2 + Qg 2 )AR 

5 

94 

-(16tt 2 ) 5Z e c = 87r(Y e *Y P T )A 2 R 2 - (—g 2 )AR 

D 

-(167T 2 ) 5Z D < = 8n(Y d *Y d T )A 2 R 2 (^g 2 + f 5 2 )Ai? 
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-(16^ 2 )<5Z Q = 4n(Y^Y u +YjY d )A 2 R 2 -(-g 2 1+ 6g 2 2 + -g 2 3 )AR 
~(l^ 2 )5 Zu , = 8n(Y:Yj)A 2 R 2 -(^g 2 + ^g 2 3 )AR 



from which we deduce the Yukawa beta functions 

1 . d 



Pv d = --A^(5Z T DC Y d + Y d 5Z Q + Y d 5Z Hd ) 
Py u = ~A-^-(6Z^ c Y u + Y u 6Z Q + Y u 6Z Hu ) 
0Y d = - l -A-^(5Z T E cY e + Y e 5Z L + Y e 5Z Hd ) 



(D2) 
(D3) 
(D4) 
(D5) 
(D6) 
(D7) 
(D8) 



2. Matter fields restricted to the brane 



There again we show all 4 diagrams contributing for the 3 generations of flavour: 



(0) 



(n) 



(0)1 



(n)l 



The sum gives 



sz™ = 



-4g 2 C 2 (R)5l 
1 I6ir 2 

*^ 2_log(Ai?) 



1 



■ log(Ai?) 8Z, 



(2) 



-8g 2 C 2 (R)S^ 
i 



16tt 2 



(4Ai?-log(Ai?)) 



sz ^ =i ^M { , AR _ log{AR)) 



SZ V ~ ~"16^ 2 



[-16Ai? 3 2 C 2 (i?)% + AAR\ lkl \* kl ] 

As for the Higgses, the gauge diagrams are those the previous case and for SZ we collect : 

1 



SZ H =- 



16tt 2 



-8ARg 2 C 2 (R) + 2 \og( A R)Tr(YiY?) 



(D9) 



(D10) 
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It is straightforward to deduce the following renormalisation constants for the matter and Higgs superfields : 



-{16tt 2 )5Z Hu = 6TT(Y^Y u )\ogKR-{-gl + Qg 2 2 )KR (Dll) 

5 

-{l^ 2 )5Z Hd = [6Tr(y d f r d ) + 2Tr(Y e ty e )] log AR- (~gl + 6g*)AR (D12) 

-(16tt 2 )6Z l = [8(YjY e )-^g 2 -12g 2 ]AR (D13) 

5 

AO 

-(16n 2 )SZ E c = [16(Y;Y C T ) -g 2 ]AR (D14) 

-(167T 2 ) sz D o = [i 6 (r;r/)-^„^ 5 2 ]Ai? (D15) 

-(l6ir 2 )6Z Q = [8(Y^Y u + YjY d )-^g 2 -12g 2 ~^g 2 }AR (D16) 

-(167T 2 ) SZuc = [16(Y:Yj)~^g 2 -^g 2 }AR (D17) 
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